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Observation of multiple ionization pathways for OCS in an intense laser field resolved 
by three dimensional covariance mapping and visualized by hierarchical ionization 

topology. 
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The two- and three-body Coulomb explosion of carbonyl sulfide (OCS) by 790 nm, 50 fs laser 
pulses focussed to ~ 10 16 Wcm -2 has been investigated by three-dimensional covariance mapping 
technique. For the first time in a triatomic molecule, a single charge state, in this case the trication, 
has been observed to dissociate into two distinct energy channels. With the aid of a three dimensional 
visualization technique to reveal the ionization hierarchy, evidence is presented for the existence of 
two sets of ionization pathways resulting from these two initial states. While one group of ions can 
be modeled using a Classical enhanced ionization model, the second group, consisting of mainly 
asymmetric channels, can not. The results provide clear evidence that an enhanced ionization 
approach must also be accompanied by an appreciation of the effects of excited ionic states and 
multi-electronic processes. 

PACS numbers: 42.50.Hz, 33.80.Gj, 33.80.Wz 



I. INTRODUCTION 

Advances in tabletop femtosecond laser generation and 
amplification techniques has lead to the routine availabil- 
ity of light pulses which may be focussed to intensities 
comparable to the binding between the proton and elec- 
tron in atomic hydrogen [lj. When a triatomic molecule 
is exposed to such a light pulse, the dynamics of electron 
removal and subsequent Coulomb explosion is rich in its 
complexity and not yet fully understood. Typically, the 
variation of the focussed laser intensity throughout the 
confocal volume 0, exposes the molecule to ~ 10 16 
Wcm~ 2 at focus, decreasing to such an extent that in 
the low intensity lobes before and after the focus that 
the effects are purely perturbative. Depending on the lo- 
cation and orientation of the molecule, the laser pulse can 
initiate a number of complex processes, including laser- 
induced reorientation low energy dissociation with 
the removal of 0, 1 or 2 electrons 0,[f| and two- or three- 
body Coulomb explosi on ( 2BCE and 3BCE respectively) 
II U H IS G3 111 13 111 El 111 El with the removal of 
at least 3 electrons. Additional complexity is introduced 
through the observation that the laser pulse may mod- 
ify the geometry of triatomics undergoing 3BCE (H2O 

pk ocs 0, co 2 0, S S El El El, so 2 E3, no 2 

|ll|). The recent development of complete momentum 
measuring coincidence techniques, E3l an d the advent 
of sub 10 femtosecond laser pulses has allowed the ob- 
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servation of the near equilibrium geometry of triatomic 
molecules through Coulomb imaging 18] . In addiction 
both concerted and sequential dissociation channels have 
been mapped in time for SOo + |Ti| ■ Concerted and se- 
quential channels in the triply charged ions of N 2 and 
CS 2 HHil have also been identified. 3BCE is fur- 
ther complicated in many-electron triatomics by the large 
number of possible ionization channels, here labelled ac- 
cording to the convention OCS( m+n+p )+ -> O m+ + C™+ 
+ S p+ and referred to as the (m, n, p) channel. 

Theoretical treatments of two- and three-body 
Coulomb explosion have been proposed and may be sep- 
arated into quantum mechanical [2! 0] and classical 
PH interpretations, broadly termed enhanced ionization 
(EI) . However recent experimental and theoretical works 
present possible shortcomings. A number of theoretical 
considerations indicate that molecular ionization cannot 
be treated as the action of a single electron. Specifically, 
ab initio calculations of the distortion of carbon dioxide 
during multiple ionization by a static electric field j2|| 
indicate the presence of a charge-exchange mechanism 
leading to experimentally observed geometries 0, and 
ionization proceeding through the most negative atomic 
site. The ongoing theoretical investigation into the appli- 
cability of Thomas- Fermi theory to ultrafast ionization in 
molecules |2jfl has lead to a recent publication concerning 
the ionization of C0 2 and N 2 0[2||. This hydrodynamic 
treatment of a group of active electrons has resulted in 
an impressive agreement between experiment and calcu- 
lation, however the extent to which the predicted dynam- 
ics are physical is open to doubt. This is particularly 
clear in the case of CI 2 [27j where the model generates 
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a CI-CI bond much stiffer than is natural which subse- 
quently plays an important part in keeping the molecule 
at low bond length as ionization proceeds. 

An important consideration in molecular ionization is 
the sequence of events within the laser pulse temporal 
envelope. A strong coupling exists between electron re- 
moval and the expansion of the molecular bonds. Thus, 
as the molecule is ionized - either sequentially or non- 
sequentially, the rate of expansion of the molecular bonds 
can vary, strongly influencing further ionization. As pro- 
posed in enhanced ionization |25j , a classical treatment of 
molecular ionization by Posthumus and coworkers, mul- 
tiple ionization is governed by dissociation, which gives 
rise to sequential ionization every time a classical ion- 
ization threshold is crossed. The 'critical separation' in 
the bond length, where the final ionization takes place, is 
simply the consequence of the strong dip in the threshold 
intensity dependence on bond length. Typically, these 
critical separations are two or three times the ground- 
state equilibrium bond length depending on the molec- 
ular type. The threshold dip is the consequence of the 
combined effect of electron localization and Stark shift- 
ing. We have previously shown that the adaptation of 
this method to triatomics such as OCS can give quanti- 
tative agreement with experiment. A criticism of both 
the classical enhanced ionization treatment and hydro- 
dynamic models is that neither makes any account the 
nonadiabatic response of a group of electrons to a laser 
pulse of a specific frequency and intensity, which has 
been shown to be influential |2^. Exact numerical so- 
lutions of the hydrogen molecular ion by Bandrauk and 
co-workers p3 | have independently predicted a maxi- 
mum in the ionization rate through a different mecha- 
nism. The electric field generated by the laser pulse pro- 
duces a Stark-shift of the lowest-unoccupied molecular 
orbital (LUMO). For particular internuclear separations, 
this Stark-shifting lifts the LUMO above the Coulomb 
barrier. However, for ionization to proceed, the LUMO 
must be populated: as the molecule expands during the 
laser pulse, the exchange of population between the high- 
est occupied molecular orbital (HOMO) and the LUMO 
becomes nonadiabatic, producing a net population in the 
LUMO. Furthermore, Kawata, Kono and Bandrauk have 
predicted similar maxima in the ionization rate in the 
case of linear H^ in an intense laser field [3(]]. As in 
the case of H J [23] , the response of the molecule to the 
laser field can be classified as either adiabatic or nonadi- 
abatic. In the adiabatic system, electron transfer along 
the molecule is possible, whereas in the nonadiabatic case 
electron migration is suppressed due to laser induced lo- 
calization. The ionization dynamics of this two-electron 
system depend heavily on both the electron distribution 
along the molecule and the nonadiabatic transitions be- 
tween the lowest three electronic states. 

These processes although not explicitly predicted for 
multi-electron systems have been the subject of some de- 
tailed experimental work |3lL l32| . By making a careful 
examination of the bond lengths at which channels of 
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FIG. 1: Schematic of the carbonyl sulfide molecule currently 
under investigation. 



the exploding N2 and I2 molecules appear to occur, it 
was deduced that certain channels could not be precur- 
sors of others and that because of this it was not possible 
that the symmetric enhanced ionization channel (1,2)- 
(2,2) could exist. It was therefore necessary to invoke 
multi-electronic processes causing asymmetric ionization 
pathways ie (1,2)-(1,3). The authors postulate that this 
process still involves the strong coupling of g — u sym- 
metry states (analogous to HOMO, LUMO in the en- 
hanced ionization symmetric ionization process) but that 
the strong coupling of the outer electron to the field gives 
rise to the ionization of an inner electron (internal re- 
scattering). Such multi-electronic processes would leave 
the molecular ion in excited states and direct evidence of 
this has been found in fluorescence studies of N2 [33L l34| 
and various other molecules pjH ]. One innovative fea- 
tures of was to introduce a visualize method able to 
show the progression of ionization through many avail- 
able pathways. 

In the case of a near-linear triatomic such as OCS, 
the molecular geometry is defined by the bond lengths 
(roc and res) and the bend angle (9) between the 
bonds, as defined in figure 1, these are the controlling 
factors in EI. The current work continues our investiga- 
tion of the behaviour of triatomics in a laser field; in the 
current work, we return to the carbonyl sulfide (OCS) 
molecule to resolve the ionization sequence. Our pre- 
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vious work |jj used momentum mapping and classical 
enhanced ionization to resolve the bond lengths at ion- 
ization for the 334 channel. We were able to show that 
the difference between the bond length calculated using a 
simplistic Coulomb assumption (the molecule is station- 
ary before explosion) is within ten percent of the correct 
bond length calculated using an iterative technique. Al- 
though ionization channels have been identified in CO2 
[3 [T3- E| , and the possibility of non-sequential ion- 
ization taking place in order to reach the (111) channel 
has been discussed [3(j the transitions between the higher 
channels have not been probed. By choosing the highly 
asymmetric OCS molecule instead of CO2 we might ex- 
pect to be able to find the limits of applicability of the 
classical ionization approach. Where possible our aim is 
to clarify the role of transitions between channels through 
the visualization of a hierarchical ionization topology. 

The covariance mapping (CM) technique was used to 
investigate laser-molecule interactions as early as 1989 
[23, but it was 1994 before exhaustive investigations into 
ionization in CO2 were published. Frasinski et al |12( em- 
ployed the technique to identify ionization channels; the 
observed channels were confirmed by Cornaggia et al [lflj 
using a variant of the technique. Cornaggia performed 
further experiments on polyatomic molecules using CM, 
reporting a straightening of carbon dioxide in 130 fs laser 
pulses flj]. Without implying any compromise in the 
conclusions published by Cornaggia and co-workers, it 
has since become apparent that the CM technique is not 
ideally suited to investigating the geometry of explod- 
ing triatomics, which requires unit collection efficiency 
for all ions, irrespective of momentum. This situation 
was addressed initially by Hishikawa and co-workers who 
proposed the mass-resolved momentum imaging (MRMI) 
technique [9I Hoi ITU ] . An analogue to this technique, ion- 
momentum imaging (IMI) , was developed by the authors 
and co-workers, and applied to the 3BCE of H 2 CO2 
and OCS where laser-induced geometry modifica- 
tions were reported and compared to Monte Carlo sim- 
ulations. Importantly, the IMI technique also includes 
an instrumental correction process , that removes 
the instrumental bias by quantifying the variation of de- 
tector efficiency with ion species and momentum. W T 
Hill and co-workers examined molecular distortion by iso- 
lating 3BCE products from a specific channel and bend 
angle, then making comparisons between the EI and the 
Thomas- Fermi models 16]. This novel approach indi- 
cates the Thomas-Fermi and EI schemes may be linked 
via an effective charge defect. 

In this paper, dissociation and 2BCE are investigated 
using two-dimensional covariance mapping (2DCM). 
Measurements of the kinetic energy of the product ions is 
reported. The 2DCM technique has previously been used 
to investigate 3BCE of C0 2 and OCS 0, however, 
given that OCS is comprised of three different atomic 
species, and that the charge-to-mass ratio of oxygen 
and sulphur atoms is degenerate under certain ionization 
states, three-dimensional covariance mapping (3DCM) of 



the Coulomb explosion of OCS is necessarily employed 
in the presented work. This allows unambiguous iden- 
tification of the ionisation channels and measurement of 
their relative strengths. The kinetic energy release asso- 
ciated with each of these channels is also reported, allow- 
ing the OC and CS bonds to be measured. This paper 
is organised as follows: the experimental configuration 
will be discussed in the following section, with details 
of the data collection and processing. Results of 2DCM 
experiments will then be presented, and details of the 
observed OCS 2BCE channels introduced. The 3DCM 
of OCS will then be presented, with particular discus- 
sion concerning the application of modern visualisation 
techniques, allowing new insights into the behaviour of 
a triatomic during Coulomb explosion. Measurements of 
the toe and r^s bond lengths will then be presented for 
each channel, along with an estimate of the relative chan- 
nel strengths. At this point, the scheme of hierarchical 
ionization topology will be introduced, and the 3DCM 
results reinter pret ed using a 3D diagrammatic technique 
adapted from |32fl . 



II. EXPERIMENTAL CONFIGURATION AND 
DATA PROCESSING 

As in previous experimental studies, the RAL (UK) 
ASTRA laser facility generated the femtosecond laser 
pulses, as detailed in previous publications H, • 
Ti: Sapphire seed pulses at 790 nm were amplified by 
CPA [38j to a pulse energy of 1 mJ in 50 fs at a repeti- 
tion rate of 10 Hz. An f/5 focus was then generated in the 
source region of a Wiley-McLaren [3^ time-of-flight mass 
spectrometer (TOFMS), interacting the laser pulses with 
the gas target. The background pressure in the TOFMS 
system is ~ 8 x 10~ 10 mbar. Product ions generated by 
the interaction of the focused laser pulse with the target 
gas were extracted from the source region, and detected 
using a pair of micro-channel plates after 11 cm of field- 
free drift. Throughout the present work, the laser polar- 
ization direction is parallel to the axis of the TOFMS. 
Unlike previous IMI experiments, the detector is now ex- 
posed to the full confocal volume and operated in high 
efficiency mode, with an extraction field of 600 Vcm -1 . 
The output of the TOFMS is monitored on a Tektronix 
TDS-744A digital storage oscilloscope (DSO). Time-of- 
flight spectra were recorded on the DSO, and stored on a 
laboratory PC for off-line analysis. The two- and three- 
dimensional covariance coefficients (C2(x, y) and Cz(x, 
y, z) respectively) were then evaluated for all n points for 
each of the N TOF spectra. For definitions of Qi(x, y) 
and Cs(x, y, z), the reader is referred to reference [40| . 
which provide full discussions of the application of co- 
variance mapping to molecular ionization. Depending on 
the process under investigation, either 02(2;, y) (2BCE) 
or 03(0;, y, z) (3BCE) was evaluated for all n points on 
the TOF spectrum, resulting in an n 2 matrix in the case 
of 02(0;, y), referred to as the two-dimensional covariance 
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FIG. 2: Two-dimensional covariance map of carbonyl sulfide 
recorded with 790 nm 50 fs Ti: Sapphire laser pulses focused 
to an intensity of ~ 10 16 Wcm -2 . The regions of positive 
covariance are identified by the labelled average time-of-flight 
spectra, the peaks are identified in table 1. The features iden- 
tified a - d are discussed in the text. 



map (2DCM) and an n 3 matrix in the case of Cs(x, y, z), 
similarly referred to as the three-dimensional covariance 
map (3DCM). Throughout this paper, the axes x, y are 
defined as the time-of-flight axes in the 2DCM, and x, y 
and z are defined as the time-of-flight axes in the 3DCM. 

To test the performance of the experimental appara- 
tus with respect to false correlations [4(| , a series of two- 
dimensional maps were recorded for the order of 3000 
laser shots as a function of target gas pressure. At a 
target gas pressure of 3 x 10 -9 mbar, the contribution 
of false correlations to the covariance map became neg- 
ligible. Of the order of 10 5 individual TOF spectra were 
recorded directly to hard disk. 



III. RESULTS - 2BCE OF OCS 

Figure 2 is a surface representation of the 2DCM of 
OCS after 10 4 laser shots, under the optical conditions 
discussed in the previous section. Correlations on the 
map are identified by labeled average time-of-flight spec- 
tra to the left and below the map, where the peaks in the 
ionization signal are identified in table I. As is apparent 
from table I, peaks 3, 5, 6, 10 and 11 are the result of the 
detection of two or more ions, either through charge de- 
generacy (peaks 5, 6, 10 and 11) or temporal overlap. As 
a consequence, the average time-of-flight spectrum gives 
a limited amount of information about the explosion of 
the molecule. Note that both TOF axes are truncated 



TABLE I: Peak assignments for the average time- of-flight 
spectra presented in the present work 



Peak Ion Species 
1 

2 0%+ 

3 C 2 + / 0%+ I S% A 

4 S|+ 

5 2 F + / S*+ 

6 O 2 / / S|+ 

7 S 3 / 

8 S 3 / 

9 C + 

10 0+ / S 2 + 

11 0+/ S 2 B + 

12 C0 F/M/B 

13 OCS 2 + 



°F/B 



between peaks 11 and 12, as no ions are detected be- 
tween t = 650 and 850 ns. The major features of figure 2 
are as follows. The strong diagonal signal is the autocor- 
relation line, a consequence of an event at x = y always 
being self-correlated and therefore gives no information 
about the processes under investigation. True correla- 
tions lie above and below this line, false correlations as- 
sociated with peak 13 (OCS 2+ ) are unavoidable. The 
double ionization of the molecular ion requires a lower 
intensity as compared to the multiple ionization leading 
to Coulomb explosion, hence the volume of the focus gen- 
erating OCS 2+ is approximately an order of magnitude 
greater than that generating the Coulomb explosion sig- 
nal 0. 

Important features of figure 2 are highlighted a - d, 
and are identified thus: (a) CO+ + S+, (b) CO+ + S 2 +, 
(c) S 3 + + 0+ and (d) S 3+ + 2 +. Correlations (a) and 
(b) identify the 2BCE channels (1, 1) and (1, 2) respec- 
tively. The kinetic energy release (KER) associated with 
these 2BCE processes, channels a and b, may be calcu- 
lated from the temporal width in x and y of the cor- 
relations, and are presented in table II. These channels, 
(1,1) and (1,2) are expected in the 2DCM, as the OC 
bond is considerably stronger than the CS bond, as ob- 
served in recent dissociative photoionization experiments 
performed by Eland and co-workers |4l| and Brion and 
co-workers 42], implying that early ionization and dis- 
sociation favour energetically efficient routes. Wang and 
Vidal 13 have published cross sections for the electron 
impact dissociative ionization of OCS, and reported the 
CO + + S + channel has a cross section at least five times 
that of any channel involving CS + . In the present work, 
only a trace of CS + is detected at a longer flight time 
than that presented in figure 2, however it is uncorre- 
cted on an equivalent 2DCM, implying + + CS + is 
negligible, and the CS + signal probably originates from 
the very weak O + CS + dissociative process. 

Because of the charge degeneracy and temporal overlap 
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FIG. 3: Geometry of the autocorrelation planes present in a 
three-dimensional covariance map 



present in much of the average timc-of-flight spectrum, 
the dense area of the 2DCM containing the atomic prod- 
uct ions where t < 0.65 (is is difficult to interpret. How- 
ever, relying on the fact that peaks 7 and 8 are purely 
due to S 3+ , correlations (c) and (d) can be employed to 
make limited inferences about the molecular behaviour. 
Correlations (c), peaks 7 (S^ + ) and 11 (O^ / S^ + ) peaks 
8 (S^ + ) and 10 (Op / S^ + ) are forwards - backwards 
pairs, formed by the ions being produced parallel to the 
laser polarization direction. Similarly in correlation (d) 
peaks 7 (S|+) and 6 (0 2 + and peak 8 (S| + ) and 5 
(0^ + /S^ + ) are forward backward pairs. As observed in 
7] the molecule favours exploding with the O - S axis 
aligned to the laser field, i.e. ip =0 in figure 1. All tri- 
atomic molecules exposed to similar laser conditions have 
been found to exhibit this forward-backward behaviour. 
However other correlations are overlapping such as those 
associated with peak 3 which represents C 2+ as well as 
3 + and S^+. 



IV. RESULTS - 3BCE OF OCS 

To establish the 3BCE channels, we turn to the 3DCM 
results. Due to the extra dimension introduced into the 
calculations, the correlation islands on the 2DCM, figure 
2, become correlation volumes in the 3DCM. True 3BCE 
correlations appear as small volumes of positive covari- 
ance within the map volume. Identifying the locations of 
these correlations is complicated by the presence of three 
autocorrelation (AC) planes, which traverse the map at 
x = y, x = z and y = z. As with the AC line in a 2DCM, 
volumes of the 3DCM located on the AC planes have a 
stronger covariance than 3BCE correlations but are es- 
sentially meaningless. The geometry of the AC planes is 
illustrated in figure 3, and will be present in any 3DCM. 

Research Systems IDL (Interactive Data Language) 
5.5 was used to examine the OCS 3DCM, allowing us 
to present, for the first time, visualizations through the 
body of the 3DCM. Previous studies [H [H Oj, HI 
present two-dimensional slices through the body of the 
3DCM, where a slice or plane parallel to the x-axis shows 
all correlations in y and z. While being a reasonable man- 
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FIG. 4: Projection visualization of the OCS three-dimensional 
covariance map. The range in (x, y, z) corresponds to the re- 
gion within the dashed box in figure 2. The peaks are identi- 
fied by comparison with table I. Three semitransparent planes 
illustrate the location of the (2, 3, 4) channel. 



ner in which to view a complex dataset, when viewed in 
three-dimensions, the location of local maxima may be 
identified more easily and accurately. A voxel (volume 
element) projection of the 3DCM between < (x, y, z) 
< 300 ns is presented in figure 4, corresponding to the 
region in figure 2 between the origin and the dashed lines. 
The full 3DCM was calculated up to t = 1 ms, but for 
the purposes of presentation, only a limited subset of this 
volume presented in figure 4. 

The visualization presented in figure 4 is generated in 
two steps. Firstly a transparency threshold, T is de- 
fined such that if C$(x, y, z) < T, the point (x, y, z) 
is treated as transparent. If 03(0;, y, z) > T, the point 
(x, y, z) is assigned a variable opacity, where the opac- 
ity depends logarithmically on G$(x, y, z). Secondly a 
three-dimensional transform is applied allowing arbitrary 
rotation, translation and oblique display of the 3DCM 
onto an image plane, in this case defined as the plane 
of the figure. Visualization in this manner is similar to 
taking an x-ray of the dataset, with the covariance coef- 
ficient equivalent to density, thus dark regions of figure 3 
correspond to volumes of high covariance. The position 
of the AC planes in figure 4 may be identified through 
a comparison with figure 3. The dark volumes between 
the AC planes are true correlations, generated only by 
3BCE channels, identified in part by average TOF spec- 
tra displayed parallel to the x, y, z axes. The labels on 
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FIG. 5: Iso-surface visualization of the OCS three- 
dimensional covariance map, where the complication of the 
dataset by the presence of the three autocorrelation planes 
has been almost negated by viewing the map along the direc- 
tion x = y — z. Again, the peaks are identified by comparison 
with table I. The circled region corresponds to the three-body 
correlation forming the (2, 3, 4) channel. 



the average TOF spectra adhere to the peak assignments 
in table I. 

An example of a 3BCE correlation between peaks 1 
(C 3+ ), 5 (0|+ / S^+) and 6 (0|+ / S*+) is identified 
in figure 4 by the intersection of three semi-transparent 
planes superimposed in the 3DCM. This correlation 
therefore represents the (2, 3, 4) channel, where the ion 
order has been changed to agree with the (m, n, p) 
convention discussed earlier. However, the KER asso- 
ciated with this channel cannot be measured, as there 
is uncertainty due to the charge-to-mass degeneracy as 
to whether axes y and z represent and S^, + respec- 
tively or vice versa: the solution of this problem will be 
presented later in this section. 

Figure 5 is an iso-surface representation of the 3DCM, 
an iso-surface being the 3D equivalent of a contour in 
2D. The image plane in figure 5 is defined normal to the 
line X = y = z, dramatically reducing the impact of the 
AC planes, as each of the AC planes x = y, x = z and 
y = z are normal to the image plane. A series of iso- 
surfaces is used to examine the 3DCM, thus identifying 
the 3BCE channels, where true correlations are identi- 
fied in the regions between the AC planes by adjusting 
the iso-surface threshold. The C3(a;, y, z) volume is in- 
vestigated by varying the surface threshold and noting 
local maxima in the regions bounded by the AC planes, 



TABLE II: Kinetic energy release (KER) associated with all 
2BCE and 3BCE channels as identified from the 2DCM and 
3D CM. The roc and res bond lengths are modeled using a 
Monte Carlo method, and the channel strength is estimated 
by integrating the covariance volumes. 

Channel KER (eV) r oc (A) r CS (A) Strength (arb.) 
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corresponding to 3BCE channels. The primary benefit 
of 3D visualization is that by moving the viewpoint from 
parallel to the x, y or z axes, an analogue to taking a 
slice, to a point on the line x = y = z simplifies the sys- 
tem considerably. Furthermore, by examining volumes of 
the map rather than slices, a more accurate identification 
may be carried out, effectively removing the AC planes 
from the dataset and thus making the volumes of interest 
immediately accessible. 

The KER associated with each of the correlations iden- 
tified is calculated from the time difference between the 
centre of the correlation volume and the zero kinetic en- 
ergy point for each ion charge-to-mass ratio parallel to 
the x, y, z axes. The total KER associated with each 
correlation volume is presented in table II, along with 
an estimate of the relative strength of the channel in the 
range to 5. The detected channel strength is gauged by 
integrating the correlation volumes. Returning to iden- 
tifying the correlations highlighted in figures 4 and 5, by 
converting the KER (eV) to momentum, p (amu ms ) 
for all combinations of ion identities possible from table 
I, a straightforward application of the conservation-of- 
momentum parallel to the detector axis allows the cor- 
rect bond length assignment to be made. An incorrect 
assignment of the 3DCM correlation volume results in a 
net momentum imbalance, whereas a correct assignment 
leads to zero net momentum. This process is repeated 
for each correlation volume in figures 4 and 5, and for all 
possible ion assignments for a particular group of three 
peaks, thus the 3BCE channels arc unambiguously iden- 
tified. The final assignments are presented in table II. 

In the case of a diatomic molecule Coulomb exploding, 
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it is trivial to estimate the bond length at the point of 
explosion 25]. However, in the case of 3BCE of OCS, 
which, as apparent from table II, is heavily predisposed 
to charge-asymmetric Coulomb explosion, such a calcula- 
tion is not straightforward. To address this, the following 
technique was applied. The bend angle (0) distribution 
between roc and res (see figure 1) for identical laser con- 
ditions is known and using the Monte-Carlo software 
developed as part of our IMI technique 0,0,0, the 3BCE 
of OCS is simulated in momentum space for a range of 
geometries. The bond lengths reported in table II are 
those in strongest agreement with the 3DCM presented 
in figures 4 and 5, with roc, r cs and 9 bound by con- 
servation of momentum. As discussed earlier these bond 
lengths assume a coulomb explosion from a stationary 
molecule. Although this is clearly not the case, previ- 
ous work suggests that even for the highest channels 
the discrepancy is around 10% between this approxima- 
tion and the correct bond length, and so rather than 
embark on the iterative process set out previously 
we will begin our analysis with these approximate val- 
ues. A number of general trends are observed in table 
II: there is an overall increase in the KER with increas- 
ing channel order, an observation common to Coulomb 
explosion studies jl H El EHH , caused by the 
higher Coulomb repulsion between the ionic constituents. 
Furthermore, the roc and rcs bonds are generally asym- 
metric throughout, and tend to increase with channel 
ionization, q = (m + n + p). Importantly, two (1, 1, 
1) channels are observed in the 3DCM, labeled (a) and 
(b) in table II in order of increasing KER. There should 
clearly be a dynamic coupling between the 3BCE chan- 
nels presented in table II, however the tabulation of the 
results is somewhat ambiguous. 



V. DISCUSSION 

To comprehend the possible pathways which the OCS 
molecule can take during ionization we first make the as- 
sumption that the ionization channels which we observe 
are the steps on an ionization ladder which finishes at 
(335). To visualize this we have adapted the method of 
[32l | into a 3D hierarchical ionization topology. The gen- 
eral hierarchical ionization (GHI) diagram for laser-OCS 
interaction is presented in figure 6, where the ionization 
channels presented in table II, identified from the 3DCM 
presented in figures 4 and 5 are indicated as nodes on a 
pyramidal framework. The least-charged (lowest order) 
channel (1, 1, 1) is the node at the apex of the pyramid, 
corresponding to a total charge, q — 3. Progression down 
the pyramid from the apex between the subsequent tri- 
angular grids indicates ionization through the loss of one 
or more electrons. Each intersection on a particular tri- 
angular grid uniquely represents a channel. Three classes 
of node are introduced in figure 6: (1) each intersection 
represents a numerically possible channel, (2) the open 
circles represent an experimentally possible channel, and 




FIG. 6: General hierarchical ionization (GHI) diagram for the 
3BCE channels observed in the 3DCM as presented in table 
II. Open circles represent numerically possible channels, solid 
circles represent true 3BCE channels. 



(3) the filled circles represent observed 3BCE channels. 
The refinement to these three classes of node makes the 
distinction between whether an ion may be generated, is 
observed in the average time-of-flight spectrum or is ob- 
served as a true 3BCE channel, following the application 
of conservation of momentum as discussed earlier. The 
key to the right of figure 6 indicates experimentally possi- 
ble channels (regular type) and observed 3BCE channels 
(bold underlined). For example, on the q = 9 level, 28 
channels arc numerically possible, illustrated by the large 
basic grid. However, the necessary component ions are 
only observed for six channels (indicated in the key), and 
of these six, only two triply-correlated channels are ob- 
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FIG. 7: The 3BCE channels identified from the 3DCM pre- 
sented in (roc, res) space, forming the second stage in hierar- 
chical ionization topology. The two groups of 3BCE channels 
observed are separated by a dashed line for clarity. 



served in the 3DCM for this level. The vertical arrows in 
figure 6 indicate the transfer between g-levels, the heads 
of which reach the vertices of three-sided secondary pyra- 
mids on the next g-level. At this point the overall charge 
of the molecule is known, but the specific ionization chan- 
nel is not. The ionization routes emanate from the vertex 
of the secondary pyramid. In choosing one of the three 
sides of this pyramid, the ionization channel is specified. 
At this point, molecular ionization is treated as a single- 
site sequential process. GHI is well illustrated by the 
transition between the (2, 2, 3) channel on the q = 7 
level and the (2, 2, 4) and (2, 3, 3) channels on the q — 
8 level. 

The most striking feature of figure 6 is the apparent 
profusion of ionization routes through the GHI diagram. 
We now turn to the second stage of hierarchical ionization 
topology, and consider the influence of the geometry of 
the exploding molecule as it develops through the GHI di- 
agram. Figure 7 shows the result of plotting in (roc, r cs) 
space all of the ionization channels uniquely identified us- 
ing 3DCM. The first observation which is clear from fig- 
ure 7 is that the approximate bond lengths are not partic- 
ularly constant as would be expected from an enhanced 
ionization explanation, although on further inspection 
they do form a group centered around extended bond 
lengths of rcs=3.7au and roc=2-7au which would cor- 
respond 2.5 and 2.25 times the equilibrium bond lengths 
which is similar to the factor of 2, observed for CO2 |^. 



The second observation is that both bond lengths in- 
crease with charge state but in an oscillating fashion (see 
shaded area figure) 0. This increasing trend is similar 
to that observed previously for CO2 [jfl. In addition the 
oscillation is accompanied by a set of channels somewhat 
offset from the others. Specifically, the group of 12 3BCE 
channels, starting with (1, 1, l)a and finishing at (3, 3, 
5) contains the majority of channels, which are generated 
at larger roc and res locations. The channels (1, 1, l)b 
to (1, 3, 4), are generated at smaller roc and res as com- 
pared to the main group. The dashed diagonal line run- 
ning from bottom left to top right serves as a visual guide, 
separating the two regions of (roc, r cs)- For illustrative 
purposes three solid lines link the ionization channels. At 
this point we must be very careful about the conclusions 
we draw. It is tempting to see the bond information, in 
table II and in figure 7, as being indicative of those bond 
lengths occupied during the steps in the ionization ladder 
illustrated in figure 6, and this is exactly the procedure 
carried out by |23, El E3 > althou gh t he bond lengths are 
calculation differently to here in |3ll l32l |. But it is im- 
portant to realize that however accurate the bond length 
measurement are, for a particular ionization channel say 
(223) , it does not mean that this ionization state uniquely 
occurs at the measured bond lengths. A higher ionization 
state, say (233), formed in a region of the focus of higher 
peak laser intensity, might display the same bond lengths 
because of enhanced ionization. In this case the molecule 
could still have passed through the (223) state, but at 
shorter bond lengths. This is best illustrated for OCS 
by figure 8 in which shows that although each chan- 
nel may ionize last at the critical distance the precursors 
to any given channel will ionize at shorter bond lengths. 
These shorter bond lengths are not measurable with a 
long laser pulse. This phenomenon is consistent with the 
recent work using short (sub 7fs) laser pulses 0, 
which clearly shows ionization at bond lengths less than 
the critical distance, if ionization is switched off rapidly 
before the critical distance is reached. For pulses of 50fs 
or longer, it is even possible that bond lengths measured 
for say the (233) channel could be shorter than those 
measured for the (223) channel but that the (223) state 
is a precursor of the (233) channel. Two 2BCE channels 
are also included in figure 7. It is assumed that the OC 
bond has not expanded during 2BCE, thus 2BCE can be 
thought of as starting from the equilibrium geometry (o, 
figure 7), then the CS bond expands while the OC bond 
is constant, undergoing ionization to the (1, 1) or (1,2) 
channels (A, figure 7). Furthermore the dissociation of 
OCS into both (1,1) and (1,2) will be quite nonCoulom- 
bic and so the approximate bond lengths given in figure 
7 are expected to be an overestimate. It is even likely 
that these channels dissociate from equilibrium. 

Figure 7 makes a distinction between two regions in 
( r oc, r cs) space. In the case of the four channels to 
the left of the dashed diagonal line, multiple ionization 
only occurs from the carbon or sulphur sites. The charge 
on the oxygen site remains 1, even though the carbon is 
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FIG. 8: Refined hierarchical ionization (RHI) diagram for the 
3BCE channels observed in the 3DCM after simplification 
as described in the text. The solid arrows indicate single 
ionization, the dashed arrows indicate double ionization. 



observed to loose up to 2 more electrons and the sulphur 
up to 3 more. In contrast, the channels to the right 
of the dashed line in figure 7 are observed to Coulomb 
explode with levels of ionization loosely bound by the 
ionization potentials of the constituent atoms. At this 
point we tentatively proposed that the (lll)b channel 
and the asymmetric channels are associated and form a 
different ionization ladder to the channels to the right of 
the line. 

Figure 8 illustrates the refined hierarchical ionization 
(RHI) diagram for the laser-OCS interaction. Visually, 
figure 8 appears simpler and better defined in comparison 
to figure 6 as there are fewer connections between nodes; 
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FIG. 9: Simulation of ionization to the right of dashed the 
line in figure 7, using the modified enhanced ionization (EI) 
model, described in the text. An analogous simulation of 
channels to the left of the line was unsuccessful. 



however, the processes implied are significant. The three 
3BCE asymmetric channels which are observed after the 
(111)6 channel are generated through two stages of dou- 
ble ionization. Dashed lines in figure 8 indicate a double 
ionization step. A pair of double ionization steps also 
occur after the (lll)a channel on route to (1,2,3) and 
(2,2,2) from (1,1,2). 

We now attempt to reproduce the results to the right 
of the dashed line in figure 7 using the results of a series 
of EI calculations. Unlike the previous calculations pre- 
sented in 0, the current calculations allow both bonds 
to expand independently, thus for each (m, n, p) channel 
an appearance intensity surface is created in (roc, r cs) 
space. This large calculation is possible by only calculat- 
ing the effective Coulomb potential (Coulombic and laser 
field) along roc and rcs, a s it has been found that this 
potential is minimized along these axes, hence EI will de- 
pend only on these conditions. This simplification only 
applies for near linear molecules (in this case 8 = 170 
degrees) parallel to the laser polarization direction, rea- 
sonable in this case given the observations of Sanderson 
et al ■ To best replicate the laser conditions responsible 
for the observations of figure 7, a wide range of laser in- 
tensities are modeled, as the laser intensity defines which 
ionization channel occurs. Each channel occupies a re- 
gion of (roc, r cs) space, indicated by the shaded areas 
in figure 9. As is apparent from comparison between fig- 
ures 7 and 9, there is reasonable agreement between the 
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measured and predicted distribution in (roc> r cs)j espe- 
cially given the number of ionization channels involved. 

As a development of the original treatment of the 
Coulomb explosion of diatomic molecules (2^, the rea- 
sonable agreement between the figures 7 and 9 is only 
possible following two major modifications to the model. 
Firstly, using the accepted ionization potentials of the 
constituent atoms, the appearance intensity minimum in 
( r oc, r cs) space (the critical geometry) is considerably 
larger than equilibrium for low q channels, and moves 
to lower rod r cs with increasing total charge q. So, 
even before taking into account the laser pulse trajec- 
tory |25j |. the trend observed in figure 7 cannot be repli- 
cated within the unmodified EI model. By decreasing 
the atomic ionization potentials with increasing charge, 
the appearance intensity minimum is varied to give a 
far better agreement with the present observations. The 
best- fit ionization potential I p ' = k I p where I p is the 
accepted ionization potential and k is the scaling param- 
eter given by k — - 0.1(m, n) + 1 or k — - 0.05p + 0.95 
depending on the channel (m, n, p). Secondly, calculat- 
ing the laser trajectory as the molecule expands reveals 
a dramatic underestimate of the bond extension for all 
channels, for example, the unmodified model predicts the 
(3,3,5) channel Coulomb explodes from roc = 2.5 A, res 
= 3.0 A. Careful investigation of the intersections be- 
tween the appearance surfaces and the laser trajectories 
over a wide intensity reveals the low q channels ionizing 
with near equilibrium bonds, thus limiting the expansion 
possible in the higher q channels. The behaviour of the 
low q channels indicates that molecular expansion is not 
purely Coulombic, a simple solution to the overestimate 
of the interaction is to reduce the charge of the ions by 
assigning fractional charges, adjusted to give a best fit to 
figure 7. It was found that the low q channels, i.e. (Ill), 
required charge suppression « 0.7, whereas the high q 
channels, i.e. (335) needed no charge suppression, inter- 
mediate channels exhibit a linear relationship between 
these limits. 

The modifications made to the enhanced ionization 
model do have valid physical basis. Firstly the modifica- 
tions to q effectively introduce a bonding aspect to the 
low charge states and make them less Coulombic. This 
deviation from pure Coulomb interaction for low charge 
states is well known in such molecules as CI2 ^3 and 
decreases with charge state. The small modification to 
the atomic IP's can be related to the ionization potential 
of the molecule, or the gap between the molecular inter- 
action potentials. For low charge states the gap is larger 
than it would be for purely Coulombic interactions be- 
cause the bonding character of the potential, represented 
by a dip below the Coulombic curve, reduces in magni- 
tude as charge state increases. With these modifications 
it is possible to represent the general features of the ma- 
jority of the ionization channels, in a way which is self 
consistent throughout and this is a significant indication 
of the overall validity of the enhanced ionization model. 
However it is equally significant that we can not repre- 



sent all final ionization channels with this approach. The 
(111)6 channel and the other channels to the left of the 
dashed line in figure 7 can not be accommodated by this 
treatment. 

In an effort to account for the remaining channels we 
should begin by explaining the two distinct energy groups 
for the (111) channel. Firstly it is sensible to consider 
enhanced ionization as it has been reasonably successful 
thus far. In the enhanced ionization picture bond expan- 
sion towards the critical distance is require before CE 
is complete. This would seem to imply the existence of 
two critical distances for the (111) channel. Although a 
secondary minimum was found in the ionization thresh- 
old curve calculated for the (334) channel previously Q 
none was found for the (111) channel. Even if a double 
structure does exist for an ionization threshold curve it 
does not easily give rise to double valued bond lengths, 
because if an ionization trajectory crosses the threshold 
curve once, the molecule will ionize and there is no ad- 
vantage to crossing more than once. A second possibil- 
ity is the presence of two distinct energetic channels in 
the break up of OCS 3+ , this might be more likely as we 
observe stable OCS 3+ in the time of light spectrum of 
this molecule, a quite unusual observation for triatomic 
molecules; it is not seen for CO2 or N2O for instance. It is 
possible that a long lived metastable trication may com- 
pletely dissociate after the laser pulse from equilibrium 
geometry. The (111)6 channel appears to have stretched 
from equilibrium but this could be deceptive just as rcs 
appears to have stretched in the (11) and (12) channel 
this may just be the effect of a non Coulombic dissocia- 
tion channel. One problem with this possibility is that 
(111)6 appears to be precursor of the other channels such 
a s (122) and this could not take place after the laser 
pulse. A third possibility is that two distinctive dissoci- 
ation channels in the OCS 2+ molecule exist, which give 
rise to different degrees of atomic excitation. The more 
highly excited molecule would ionize to the (111) channel 
at a lower laser intensity, this phenomenon has been ob- 
served for metastable argon ions [iEj. A consequence of 
this would be to shift the critical distance to lower bond 
length, in fact two curves would exist for the two differ- 
ent molecular states with the same charge. It should be 
notes that two distinct kinetic energy groups have been 
observed for the S fragment in the CO+S dissociation 
of OCS after 233nm absorption |4^. This phenomenon 
could not be predicted from the simple classical mode, 
as it would require two distinct sets of atomic ionization 
potentials to exist. Excited molecules formed in high in- 
tensity regions would reach the (111) channel earlier in 
the laser pulse and this may be why asymmetric double 
ionization is possible after the (111)6 channel is formed, 
as multi electronic asymmetric processes occur at bond 
lengths closer to equilibrium than the critical distance 
32]. If asymmetric ionization is preferred at shorter bond 
lengths then only the state which gives rise to the (111)6 
channel can make the triply ionized molecule at suffi- 
ciently small bond length to allow further asymmetric 
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multi-electronic processes to occur. The state which give 
rise to the (lll)a channel does not do so until after the 
molecule has passed through this bond length range. 

VI. CONCLUSION 

The kinetic energy release during two- and three-body 
Coulomb explosion of carbonyl sulfide has been mea- 
sured using three-dimensional covariance mapping. For 
the first time in a triatomic molecule two distinct en- 
ergy groups have been identified for the same dissociation 
channel. A 3D hierarchical ionization topology has been 
introduced to help visualize the multitude of possible ion- 
ization pathways and a Classical trajectory simulation 
using a simple over the barrier enhanced ionization model 
has been used in an attempt to reproduce the observed 
bond lengths. The distribution of most bond lengths are 
reasonably predicted when small modification to the en- 
hanced ionization model are introduced. The presence of 
the two (111) channels and highly asymmetric channels 



however can not be modeled using the simple approach. 

In conclusion, OCS has proven to be a useful subject 
of study in the quest to understand the ionization mech- 
anisms of small molecules in an intense laser pulse. The 
results indicate the need to develop an enhanced ioniza- 
tion model which incorporates molecular excitation and 
the possibility of multi-electronic processes in order to 
be able to fully predict the ionization process. In addi- 
tion the need for further experimental study is indicated, 
in order to better understand the interplay between the 
ionization routes indicated by the current work, utiliza- 
tion of sub 10 femtosecond pulses would enable prob- 
ing within the critical distance and full coincidence tech- 
niques would establish the link between geometry and 
ionization pathway. 
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